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Background. Histological behavior of glioblastoma multiforme suggests it would benefit more from a global rather 
than regional evaluation. A global (whole-brain) calculation of diffusion tensor imaging (DTI) derived tensor metrics 
offers a valid method to detect the integrity of white matter structures without missing infiltrated brain areas not seen 
in conventional sequences. In this study we calculated a predictive model of brain infiltration in patients with glioblas- 
toma using global tensor metrics. 

Methods. Retrospective, case and control study; 1 1 global DTI-derived tensor metrics were calculated in 27 patients 
with glioblastoma multiforme and 34 controls: mean diffusivity, fractional anisotropy, pure isotropic diffusion, pure 
anisotropic diffusion, the total magnitude of the diffusion tensor, linear tensor, planar tensor, spherical tensor, relative 
anisotropy, axial diffusivity and radial diffusivity. The multivariate discriminant analysis of these variables (including age) 
with a diagnostic test evaluation was performed. 

Results. The simultaneous analysis of 732 measures from 12 continuous variables in 61 subjects revealed one discrimi- 
nant model that significantly differentiated normal brains and brains with glioblastoma: Wilks' A = 0.324, x 2 (3) = 38.907, 
p < .001 . The overall predictive accuracy was 92.7%. 

Conclusions. We present a phase II study introducing a novel global approach using DTI-derived biomarkers of brain 
impairment. The final predictive model selected only three metrics: axial diffusivity, spherical tensor and linear tensor. 
These metrics might be clinically applied for diagnosis, follow-up, and the study of other neurological diseases. 
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value of tests 



Introduction 

Some pathologic and magnetic resonance (MR) im- 
aging characteristics of astrocytomas grades II to IV 
(highest degree known as glioblastoma multiforme, 
GBM) suggest these tumors would benefit from the 
use of a global measurement of brain impairment. 1 
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The first imaging approaches to characterize 
high-grade glial lesions, especially GBM, were 
fraught with pitfalls resulting from the marked 
heterogeneity of both glial-infiltrated and normal 
brains. 23 These tumors frequently contain multiple 
areas of variable histologic features, so that a sam- 
pling error in a biopsy may mean that the degree of 
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malignancy seen by the neuropathologist may not 
reflect the degree of malignancy present elsewhere 
in the tumor, resulting in significant undergrading 
of some lesions. 2 Thus, even when all radiologi- 
cally visible portions of a tumor have been excised, 
the surgical margins may not be "clean", and fur- 
ther neoplastic growth can (and usually does) oc- 
cur in the adjacent brain tissue, leading from mi- 
croscopic residual to gross recurrence. 4 Therefore, 
none of the MR protocols for GBM in every day 
practice should be only morphologic. 3 ' 5 6 As a con- 
sequence, surgery usually only reduces the tumor; 
this information is relevant as recent evidence has 
proved gross total resection (surgical margin sta- 
tus) significantly correlates with progression-free, 
recurrence pattern and overall survival in patients 
with GBM. 7 8 

GBM is considered a whole brain disease. 
Radiotherapy and chemotherapy follow surgery. 

Although MR perfusion and spectroscopy 9 , 
and sometimes diffusion tensor imaging (DTI) 10 
are routinely used methods to locate parts of the 
tumor-GBM with high malignancy that should be 
biopsied, the development of specific and sensitive 
biomarkers remains a critical unmet need. 11 

Our purpose in this study was to explore the 
diagnostic ability of a global (whole brain) assess- 
ment of DTI-derived tensor metrics in normal and 
infiltrated brains with GBM. We used the multivar- 
iate technique of linear discriminant analysis (DA), 
previously reported in MRI diagnosis 12 , to classify 
the study participants into groups, describe group 
differences and to assess the relative importance of 
DTI variables for discriminating between groups. 
This analysis might unveil findings and associa- 
tions that cannot, in a partial-regional assessment, 
be recognized at surgery, neurologic, MRI and/or 
pathologic examination. Considering there is still 
scarce information in the medical literature about 
the global calculation of tensor metrics 1314 , a pre- 
dictive discriminating model may offer an inno- 
vative diagnostic approach to the surgical-neuro- 
oncology team. 

Subjects and methods 

Subjects 

This was a case-control study. We included pa- 
tients with suspected diagnosis and later patho- 
logical confirmation of primary GBM who had un- 
dergone preoperative brain MR examinations be- 
tween January 2010 and September 2012. Exclusion 
criteria were corticosteroid or antibiotic treatment, 



lesions with areas related to calcification and/or 
hemorrhage and previous brain surgery. A control 
group included young and elderly healthy volun- 
teers recruited from the enrolled interns and medi- 
cal residents of the hospital, and elderly subjects 
from our Geriatrics unit. All volunteers received 
detailed health examinations; exclusion criteria 
were major neurological, psychiatric, or cardiovas- 
cular diseases. A radiologist interpreted the MR 
images blinded to the patient's history and MRI 
examinations with other structural abnormali- 
ties were excluded. The local institutional review 
board approved the study (Project #2011-EXT-05). 

Brain image acquisition 

MR sequences included conventional axial 
T2-weighted imaging, axial Fluid-Attenuated 
Inversion Recovery (FLAIR), axial Spoiled 
Gradient Echo (SPGR), DWI and axial Tl-weighted 
imaging, using 0.1 mmol/kg of body weight of ga- 
dopentetate dimeglumine (Magnevist; Schering, 
Berlin, Germany). Healthy volunteers did not re- 
ceive endogenous contrast. DTI was performed 
using a single-shot SE EPI sequence. Diffusion 
gradients were applied in 25 directions with b-val- 
ues of 1000 s/mm 2 and an image without diffusion 
weighting with b-value of 0 s/mm 2 . DTI sequences 
were acquired in the axial plane with 44 contiguous 
sections, 2.4 mm section thickness, no intersection 
gap; TR/TE of 17,000/80 ms, with parallel imaging 
to reduce off-resonance artifacts (PI factor was 2); 
25 x 25 cm FOV; and 128 x 128 matrix/pixel size. 
MR was performed on a single occasion using a 3T 
unit (Signa HDxt, GE Healthcare, Waukesha, WI, 
USA); and a high-resolution eight-channel head 
coil (Invivo, Gainesville, FL, USA). 

Image postprocessing and data analysis 

We used the software dcm2nii 15 and the FMRIB 
Software Library (FSL) v. 4.1.9. 16 DTI images were 
extracted using the Brain Extraction Tool (BET) 
v. 2.1. 17 Eddy currents were corrected using the 
FMRIB 's Diffusion Toolbox v. 2.0; the Reconstruct 
Diffusion Tensor (DTIFIT) and the fslmaths tool gener- 
ated the eigenvector and eigenvalue maps for each 
tensor metric. The fslstats tool calculated the scalar 
measures (mean values) of each whole-brain cal- 
culation. The apparent diffusion coefficient (ADC) 
value, a simple index calculated from diffusion- 
weighted images 18 , was considered equivalent to 
the MD (mean diffusivity) metric, as it was ob- 
tained from the DTI sequence. 19 DTI-derived tensor 
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TABLE 1. Correlations of tensor metrics, controlled for the effect of diagnosis, age and gender 



Tensor metric 



Cs 



FA 


Pearson's 
R 

p-value 


-.552 
< .001 


FA 














RA 


Pearson's 
R 

p-value 


-1.000 
< .001 


.557 
< .001 


RA 












Cp 


Pearson's 
R 

p-value 


-.937 
< .001 


.584 
< .001 


.912 
< .001 


Cp 










CI 


Pearson's 
R 

p-value 


-.898 
< .001 


.541 
< .001 


.943 
< .001 


.673 
< .001 


CI 








L 


Pearson's 
R 


.21 1 


.075 


-.175 


-.191 


< .001 


L 








p-value 


.165 


.596 


.256 


.213 


.999 










Pearson's 


.195 


-.079 


-.183 


-.194 


.006 


.890 






P 


R 














P 






p-value 


.205 


.580 


.240 


.214 


.972 


< .001 






AD 


Pearson's 
R 


.034 


.209 


-.008 


-.002 


.106 


.882 


.880 


AD 




p-value 


.826 


.137 


.958 


.989 


.508 


< .001 


< .001 




MD 


Pearson's 
R 


.195 


-.078 


-.183 


-.193 


.007 


.892 


1.000 


.881 

MD 




p-value 


.206 


.589 


.241 


.214 


.968 


< .001 


< .001 


< .001 


RD 


Pearson's 
R 


.213 


-.226 


-.209 


-.149 


-.105 


.815 


.973 


.779 .973 

RD 




p-value 


.151 


.103 


.163 


.316 


.502 


< .001 


< .001 


< .001 < .001 


q 


Pearson's 
R 


-.306 


.804 


.310 


.339 


.328 


.403 


.281 


.627 .284 .214 




p-value 


.031 


< .001 


.030 


.016 


.026 


.003 


.046 


< .001 .044 .116 


AD = axial diffusivity; CI = 


linear tensor; Cp 


= planar tensor; Cs = 


spherical tensor; FA 


= fractional anisotropy; L 


= the total magnitude of the diffusion tensor; MD = me 



p = pure isotropic diffusion; q = pure anisotropic diffusion; RA = relative anisotropy; RD = radial diffusivity 



metrics formulas using the major (Al), intermediate 
(A2), and minor (A3) eigenvalues allowed the cal- 
culation of the eleven most common tensor metrics 
for brain imaging: mean diffusivity (MD), fractional 
anisotropy (FA), pure isotropic diffusion (p), pure 
anisotropic diffusion (q), the total magnitude of the 
diffusion tensor (L), linear tensor (CI), planar tensor 
(Cp), spherical tensor (Cs), relative anisotropy (RA), 
axial diffusivity (AD) and radial diffusivity (RD) 10 ; 
each one representing a single global measure of 
the whole-brain. Figure 1 A shows the algorithm for 
measuring the DTI-derived tensor metrics. 

Statistical analysis 

Study design 

The study was considered a Phase II aimed to de- 
termine the capacity of DTI-derived biomarkers to 



distinguish between people with cancer and those 
without. 20 

Sample size 

Considering our predictive model to discriminate 
between normal brains vs. brains infiltrated with 
GBM underwent a diagnostic performance assess- 
ment, the adequacy of the sample size to expect 
validity from our results was based on matching 
this phase with the summarized list of computed 
sample sizes needed for an exploratory retrospec- 
tive study reported by Obuchowski et al. 11 , at least 
10 diseased patients and 10 control patients were 
required to maintain statistical validation in a di- 
agnostic test evaluation where the type I error rate 
was set at 0.05, type II error rate was < 0.10, and 
power > 0.90. Our study included 27 patients and 
34 controls. 
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FIGURE 1. (A), FSL software algorithm used in the image postprocessing and data analyses. (B-E), Examples of acquired sequences in a patient with 
GBM and the tensor-metric maps generated for the data analyses: (B), axial T2-weighted; (C), post contrast axial Tl -weighted; (D), axial diffusivity (AD) 
tensor map; and (E), fractional anisotropy (FA) tensor map. Notice how it might not be possible to perform an imaging diagnosis based only on a visual 
inspection of these maps. 



Multivariate DA 

We ran a DA, which was optimal under the same 
conditions where Manova was optimal; then at- 
tempted to detect any deviation from Manova as- 
sumptions that might distort the tests of statistical 
significance. 22 We assessed the normality of the 
distribution of the DTI-derived scores using the 
Kolmogorov-Smirnov's and Shapiro- Wilk normal- 
ity tests 23 ; eliminated significant outliers, evaluated 
multivariate normality and linearity, and tested 
the homogeneity of variance-covariance matrices 
using the Box's M test. 24 Considering the similarity 
of the tensor-metric formulae, we ran scatterplots 
and correlations to check the strength of correla- 
tions among the dependent variables in order to 
detect the presence of multicollinearity and singu- 
larity (Table 1). Partial correlation analyses were 
carried out to calculate the Pearson's correlation 
coefficient (r) controlling for the effect of age, gen- 
der and clinical diagnosis. The strength of the lin- 
ear relationship corresponding to each correlation 
coefficient value was interpreted as very strong (at 
least of 0.8), moderately strong (0.6 up to 0.8), fair 
(0.3 up to 0.6) and poor (less than 0.3). A squared r 
value represented the coefficient of determination, 
the proportion of variance that each two compared 
variables had in common. 25 



We applied the stepwise method in DA, it con- 
sidered the value of Wilk's lambda and changing 
criteria: minimum partial F to enter of 3.84 and 
minimum partial F to remove of 2.71. 22 Continuous 
variables were included with the predictive aim to 
identify specific tensor-metric attributes in GBM 
and normal brains. The dependent variable (DV) 
used in the DA was the clinical diagnosis, which 
classified subjects as patients or controls. The in- 
dependent variables (IVs) included 11 DTI-derived 
tensor metrics: MD, FA, p, q, L, CI, Cp, Cs, RA, RD 
and AD, and the patients' age (in years). The effect- 
size measure for discriminant analysis was calcu- 
lated using the squared canonical correlation as the 
equivalent of the R 2 in regression. 26 By convention, 
effect sizes of 0.02, 0.15, and 0.35 are termed small, 
medium, and large, respectively. 27 For all analyses, 
statistical significance was indicated by a p-value < 
0.05. 

Diagnostic model evaluation 

The cross-validated contingency Table generated 
by the DA was used to evaluate the diagnostic per- 
formance of the DA model. We reported values of 
sensitivity, specificity, positive and negative likeli- 
hood ratios, and positive and negative predictive 
values, with their corresponding confidence in- 
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TABLE 2. Multivariate analysis (between-groups) of diffusion tensor imaging (DTI) -derived tensor metrics and age showing the statistical differences 
between means of normal-brain and brain-with- as glioblastoma multiforme (GBM) groups for the independent variables included in the analysis 



Variable 


Healthy brains 
Mean SD 


Brains with GBM 
Mean SD 


Wilks' 

LUiIlUUU 


Ftest 


p-value 


Cs (spherical tensor) 


.747091 


.026938 


.768395 


.042299 


.915 


3.341 


.076 


FA (fractional anisotropy) 


.287029 


.011517 


.254082 


.026761 


.607 


23.341 


< .001 


RA (relative anisotropy) 


.233436 


.025491 


.209370 


.033559 


.855 


6.088 


.018 


Cp (planar tensor) 


.138366 


.013823 


.136635 


.036218 


.999 


.036 


.850 


CI (linear tensor) 


.114543 


.013991 


.098463 


.011930 


.711 


14.621 


.001 


L (total magnitude of the diffusion tensor) 


.002277 


.000087 


.002117 


.000147 


.691 


16.077 


< .001 


p (pure isotropic diffusion) 


.002107 


.000077 


.001959 


.000134 


.681 


1 6.893 


< .001 


AD (axial diffusivity) 


.001548 


.000044 


.001399 


.000087 


.461 


42.052 


< .001 


MD (mean diffusivity) 


.001217 


.000044 


.001132 


.000078 


.685 


16.526 


< .001 


RD (radial diffusivity) 


.001051 


.000050 


.000997 


.000078 


.852 


6.237 


.017 


q (pure anisotropic diffusion) 


.000452 


.000036 


.000367 


.000047 


.483 


38.529 


< .001 


Age 


40.333 


21.502 


47.150 


15.187 


.965 


1.294 


.263 



SD = standard deviation 



tervals (CI). Evaluation of the diagnostic tests fol- 
lowed the Standards for Reporting of Diagnostic 
Accuracy (STARD) initiative. 28 

Software 

All analyses were carried out using the IBM® 
SPSS® Statistics software (version 22.0.0.0 IBM 
Corporation; Armonk, NY, USA). Diagnostic per- 
formance was assessed using MedCalc® (version 
12.3.0.0 MedCalc Software bvba, Mariakerke, 
Belgium). 

Results 

Subjects and MRI acquisition 

The study was conducted in 61 subjects; 27 pa- 
tients: 13 females (mean age 50.0 ± 15.400 years, 
range 31-73 years) and 14 males (mean age 46.93 ± 
15.403 years, range 18-78 years); and 34 controls: 26 
females (mean age 41.04 ± 22.37 years, range 21-80 
years) and 8 males (mean age 42.88 ± 21.89 years, 
range 24-72 years). The eleven DTI tensor-maps 
plus the age (per subject) added up 732 measure- 
ments included in the analyses. Figure 1 B-E shows 
examples of some of the MR sequences and tensor- 
metric maps used in the data analyses. 

Partial correlation analyses 

A scatterplot showed no serious violation of the 
assumptions of linearity, homoscedasticity, and 



outliers. Among 55 pairs of bivariate correlations, 
we found only 15 with a significantly very strong 
(at least 0.8) r value: Cs^RA (-), Cs^Cp (-), Cs^L 
(-), FA^q (+), RA^Cp (+), RA^Cl (+), L^p (+), 
L^AD (+), L^MD (+), L^>RD (+), p^AD (+), 
p^MD (+), p^RD (+), AD<=>MD (+), and MD^RD 
(+). Table 1 and Figure 2 depict correlation values 
and the scatterplot of the eleven tensor-metrics. 

Discriminant analysis 

Although some r values were calculated at > 0.8, 
we included all variables in the DA, as we found 
evidence the stepwise variant of this method pro- 
tects against multicollinearity and singularity; 29 a 
brief explanation is presented in the discussion sec- 
tion. The assumption of homogeneity of variance- 
covariance matrices was interpreted as significant 
(Box's M value = 35.110, F = 5.317, df (6, 9087.738), 
p = < .001). In the stepwise statistics, at each step, 
the best variable that minimized the overall Wilks 7 
Lambda was entered: AD was entered at Step 1, F 
= 42.052 (1, 36) p < .001; CI entered at Step 2, F = 
29.609 (2, 35) p < .001; and Cs entered at Step 3, F = 
23.672 (3, 34) p < .001. 

DA revealed one discriminant function that 
significantly differentiated the normal brains and 
GBM brains: Wilks 7 A = 0.324, x2 (3) = 38.907, p < 
.001. By indicating the significance of the discrimi- 
nant function, Wilks 7 lambda provided a moderate 
proportion of total variability not explained by the 
model of 10.49%. A canonical correlation of 0.822 
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FIGURE 2. Scatter matrix of the data variables grouped by diagnosis. 



suggested the model explains 67.56% of the varia- 
tion in the grouping variable. 

Summary of discriminant functions 

The tests of equality of group means provided sta- 
tistical evidence of significant differences between 
means of normal brains and brains with GBM in 9 
of the IVs, with AD producing the highest F's val- 
ue; Table 2 depicts the means, standard deviations 
(SD) and F's tests values (between-groups multi- 
variate analysis). 

Standardized canonical discriminant function coef- 
ficients showed an index of the importance of each 
predictor for diagnosis with the sign indicating 
the direction of the relationship. A significant in- 
crease in values of Cs (spherical tensor), CI (linear 
tensor) and AD (axial diffusivity) were the strong- 
est diagnostic predictors. The variable coefficients 
stood out (for these data) as those that strongly 
predicted allocation to the normal-brain or tumor- 
brain group. The coefficient score decrement was 
proportional to less successful diagnostic predic- 
tors (Table 3A). 

Structure Matrix Data provided another way of 
indicating the relative importance of the diagnos- 
tic predictors by showing the correlations (Pearson 



coefficients) of each variable with each discrimi- 
nate function. Many researchers consider the 
structure matrix correlations more accurate than 
the standardized canonical discriminant function 
coefficients. 26 By identifying the largest loadings 
for each discriminate function, different patterns of 
loading variables can be seen. We found AD, MD, 
p, q, CI, RD and FA, as the functions that best dis- 
criminate between normal brains and brains with 
tumor. A value of 0.30 was considered as the cut- 
off between important and less important variables 
(Table 3B). 30 

The canonical discriminant function coefficients 
Table showed the unstandardized coefficients (b) 
that were used to create the discriminant func- 
tion (equation), they operated just like a regression 
equation, allowing us to build a predictive model 
of brain status: 

Brain status (normal brain vs. tumor infiltration) = 
-48.295 +11443.557 (axial diffusivity, AD) +105.124 
(longitudinal tensor, CI) +26.804 (spherical tensor, Cs) 

The discriminant function coefficients (b) indi- 
cated the partial contribution of each variable to 
the discriminate function controlling all other vari- 
ables in the equation (Table 3C). 

The group centroids values described each group 
in terms of its profile, using the group means of 
the predictor variables called centroids. The cut- 
off value was defined as the mean of the two cen- 
troids; if the discriminant score of the function of a 
new case was less than or equal to the cut-off, the 
case was classed as 1 (brain with tumor), whereas 
if it was above the cut-off, it was classed as 0 (nor- 
mal brain). In our study, normal brains had a mean 
of 1.483 while brains with GBM produced a mean 
of -1.334; the cut-off for the function at group cen- 
troids showed a calculated value of 0.149. 

For the final part of the DA we performed a clas- 
sification phase using the cross- validated set of data 
to present the power of the discriminant function. 
These results revealed that 92.7% of patients were 
classified correctly into "normal brain" or "brain 
with GBM" groups, this value corresponded to 
the overall predictive accuracy of the discriminant 
function. Additional results of diagnostic tests per- 
formance including the 95% confidence intervals 
(C.I.) showed: sensitivity = 100.00 (80.49 - 100.00); 
specificity = 87.50 (67.64 - 97.34); (+) likelihood ra- 
tio = 8.00 (82.78 - 23.06); (-) likelihood ratio = 0.00 
( - ); (+) predictive value = 85.00 (62.11 - 96.79); and 
(-) predictive value = 100.00 (83.89 - 100.00). 

The average discriminant (D) scores for each 
group and the group centroids were used as vis- 
ual demonstrations of the effectiveness of the dis- 
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TABLE 3. Independent variables included in the discriminant analysis. A, ordered by their Standardized Canonical Discriminant Function Coefficients 
(variables with larger coefficients stand out as those that strongly predict allocation to each diagnosis). B, Within-groups correlation matrix depicts 
the participant variables ordered by absolute size of correlation (Pearson coefficients) within function. The largest loadings for each discriminate 
function (AD was the largest) suggest the preference of diffusivity values that discriminates between normal- and brain-tumor groups. A value of 0.30 is 
considered as the cut-off between important and less important variables, notice that variables with (*) were not used in the analysis. C, unstandardized 
coefficients used to create a discriminant function operating just like a regression equation. Coefficients indicate the partial contribution of each 
variable to the discriminate function controlling for all other variables in the equation 



A 


B 




C 


Standardized Canonical Discriminant 
Function Coefficients 


Structure Matrix 




Canonical Discriminant Function 
Coefficients 


Function 


Variable 


Function 


Function 


Variable 

1 


1 


Variable 

1 



CI (linear tensor) 


1.361 


AD (axial diffusivity) 


.748 


AD (axial diffusivity) 


11443.557 


Cs (spherical tensor) 


.962 


MD (mean diffusivity)* 


.568 


CI (linear tensor) 


105.124 


AD (axial diffusivity) 


.806 


p (pure isotropic diffusion)* 


.566 


Cs (spherical tensor) 


26.804 






L (total magnitude of the diffusion tensor)* 


.553 


(Constant) 


- 48.295 






q (pure anisotropic diffusion)* 


.533 










CI (linear tensor) 


.441 










RD (radial diffusivity)* 


.427 










FA (fractional anisotropy) 


.320 










RA (relative anisotropy)* 


.278 










Cs (spherical tensor) 


-.211 










Age* 


.124 










Cp (planar tensor)* 


.020 







criminant function. Histograms and box plots of 
the average D scores for each group were used as 
graphical demonstrations of the effectiveness of 
the discriminant function, the absence of overlap 
of the plots revealed an excellent discrimination 
(Figure 3 A-B). 

Discussion 

The lack of consensus regarding which DTI- 
derived tensor metrics are the most meaningful 31 , 
and the scarce information about their diagnostic 
abilities, compelled us to evaluate whether a global 
approach might have clinical applicability. We con- 
sider our study an introduction to the method and 
proof-of-principle that a global approach using se- 
lected DTI-derived tensor metrics can differentiate 
normal brains from brains infiltrated with GBM, 
the selected metrics may function as biomarkers 
assembling a predictive model of tumor infiltra- 
tion. 

The relevant findings in our study showed that a 
multivariate DA of global measurements excluded 
the pair-wise comparisons from conventional tu- 
mor-region evaluations; the assembled statistically 



significant discriminant model of tumor brain im- 
pairment (for these data) needed only three global 
DTI-derived metrics: AD, CI, and Cs. 

Some advantages of a global approach using DTI 
metrics need to be mentioned: it decreases the bias 
associated with manual placement of a region of in- 
terest encompassing tumor regions; the tumor and 
edema regions are implicitly included in the evalu- 
ation; lesions not perceived by the radiologist's eye 
on conventional sequences would be included in 
a global assessment; it may avoid problems asso- 
ciated with partial volume effects, and inaccurate 
image coregistrations; DTI biomarkers can be ap- 
plied to other tumors/neurological diseases; its ac- 
quisition does not need contrast, and its post pro- 
cessing method can be semiautomatic; these facts 
broad the clinical applicability with no significant 
increase in the cost of MRI examinations. 

The selected biomarkers in our final model de- 
serve a brief explanation: AD depicted the main in- 
fluence (larger value of its b unstandardized coef- 
ficient); it represents the directional diffusivity de- 
scribing the microscopic water movement parallel 
to axonal tracts. AD is one of the best biomarkers in 
the diagnosis of enhancing rim in GBM, but not for 
other tumor regions 10 ; (this fact provides evidence 
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that a regional measurements may not be the most 
effective way to use DTI metrics in brain tumor im- 
aging). 14 AD has been studied in animal models of 
encephalomyelitis of the spinal cord 32 33 , in unfixed 
ex vivo human brains with multiple sclerosis 34 , in a 
model of axonal injury caused by stroke 35 , and in 
optic neuritis. CI and Cs on the other hand, along 
with FA, have been reported among the biomark- 
ers with best overall performance in differentiat- 
ing the cystic cavity in abscess from GBM. 13 They 
show best diagnostic performance in the detection 
of normal-appearance white matter (NAWM) and 
the cystic cavity in brains with GBM. 10 Cs repre- 
sents spherical normalized coordinates of a non- 
orthogonal DTI-derived tensor for each voxel, and 
CI corresponds to the linear case. 36 
Several limitations in this study need to be ad- 
dressed: because there have not been studies inves- 
tigating a whole set of tensor metrics (not only FA 
and ADC) in a global approach 1 ' 3738 , it is difficult to 
compare our results with others in the literature. 
Further studies might include comparisons with 
other brain tumors, the influence of variables like 
radiation necrosis, inflammatory and demyelinat- 
ing diseases; and tumor infiltration categories such 
as post-surgery and post-radiotherapy; all of them 
were beyond the scope of this study. A concern 
of using DTI-metric values with high correlations 
(correlations up 0.8 or 0.9), as we observed in our 
data, might be raised because in those situations 
one variable is a near-linear combination of the 
other variable (the variable provides information 
that is redundant to the information available in 
one or more of the others, making matrix inversion 
unreliable). 29 The usual solution is a deletion of the 
redundant variable, however, because we have a 
compelling theoretical reason to retain all variables 
in this study (to evaluate the simultaneous discri- 
minant ability of 11 global tensor metrics), the IBM® 
SPSS® Statistics software protects against multicol- 
linearity and singularity through computation of 
pooled within-cell tolerance (1-squared multiple 
correlation, SMC) for each variable. SMC is the 
squared multiple correlation of a variable where it 
serves as the dependent variable (DV) with the rest 



FIGURE 3. Visual demonstration of the effectiveness of the 
discriminant function. (A), histograms showing the distribution 
of discriminant scores for normal- and tumor-brains. (B), box 
plots of the average D scores. Both kinds of plots illustrate the 
distribution of the discriminant function scores for each group. 
The box-plots depict a visual demonstration of the excellent 
discrimination of the model by showing no overlap between 
groups. 



Radiol Oncol 2014; 48(2): 127-136. 



Roldan-Valadez E et al. / Global diffusion tensor imaging in glioblastoma multiforme 



135 



as independent variables (IV) in multiple correla- 
tion. Variables with insufficient tolerance are delet- 
ed from the analysis; this procedure is a part of the 
stepwise method in DA. 22 Our discriminant model 
was able to explain a significant proportion of the 
variability in the data (67.56%), but may still have 
some errors in predicting individual diagnosis, so 
model validation should be done in subsequent 
studies. We acknowledge the linking of tensor- 
metric values with the axonal-integrity status rep- 
resents an oversimplification with respect to what 
is happening in brains with GBM, where complex 
tissue changes occur and affect water diffusivity: 
density of fiber, average diameters, degree of my- 
elination, directional similarity, cellularity, viscos- 
ity, permeability, and histologic architecture; the 
DTI-tensor values are the effects of the summation 
of all these microstructural barriers. 13 

Several questions remain unanswered, for ex- 
ample, what is the relation of these tumor-DTI bio- 
markers with those of MR perfusion and spectros- 
copy? What is the association of DTI-biomarkers 
with the pattern of relapse and extension of resec- 
tion in GBM? So far, only one study, to the best 
of our knowledge, has correlated a few regional 
DTI-tensor metrics with the survival of patients 
with GBM 39 ; thus the clinical value of global DTI- 
metrics in predicting the overall survival has yet to 
be determined. As a phase II study, our research 
line will look for a sequel, applying the proven 
concepts in the follow-up of tumor-infiltration cat- 
egories (post-surgery, post-radiotherapy, etc.) and 
in differential diagnoses (primary brain tumors vs. 
metastasis vs. demyelinating diseases). 

Conclusions 

Although we cannot affirm the superiority of glob- 
al vs. regional DTI-derived tensor metrics in the 
evaluation of GBM yet, we can ascertain with cer- 
tainty that there is an immediate clinical applicabil- 
ity of these biomarkers in assembling statistically 
significant predictive models able to announce the 
conversion of normal tissue to tumor infiltrated tis- 
sue before the conventional MR sequences show 
conspicuous findings. These principles could eas- 
ily be extended to other neurological diseases. A 
first step in the advanced evaluation of brain tu- 
mors might include a global measurement of DTI- 
biomarkers able to pick up major infiltration zones. 
Due to the large number of variables (qualitative 
and quantitative) that must be analyzed in con- 
temporary brain MRI by radiologists and neuro- 



scientists conducting research on novel imaging 
biomarkers; multivariate techniques, like DA, may 
help in the generalization of knowledge beyond 
one setting. 
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